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ABSTRACT 


The accuracy of proposed lower confidence limits for system 
availability is analyzed. Random values of the lower 100(1- « )% 
eonmdence limit Ag ee) for system availability are computed for 
a system whose failure density is exponential ( 1< and whose repair 
density is exponential VAL ). The system is modeled as an alternating 
renewal process. The 1O0O0(l-x )th percentile point of the generated 
distribution of Re iy ) is compared with system availability as a 


measure of accuracy for Ag Li)’ 
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loeiINLO@DU CTION 


The system considered in this thesis is of the type which operate 
for a time, fail, are repaired and returned to operation. The system 
then may be said to have two states and can be modeled as a renewal 
peocess. 

The objectives of this paper were the following: 

1. To compute exact values of the steady state availability for 
a system with given distributions of times to failure and times to 
repair. 

2. To compare the availability of a system based on the test 
procedures used with the 100(1-¢) th percentile point of the dis- 
tribution of a proposed 100(1 - « )% lower confidence limit. This 
will allow an accuracy analysis of a proposed lower confidence limit. 

The method used was to fit a normal distribution to simulated 
data and compute the 100(1- ) th percentile point of the distribution. 

An alternative method was investigated, which fitted a 2 para- 


meter gamma distribution. The procedure proved to be inaccurate. 





Il STATEMENT OF THE PROBLEM AND ANALYSIS 


A, DEFINITION OF AVAILABILITY 

In this paper the term availability [1] is defined as the probability 
that a system will be "up" or operable when called upon at a random 
time. 

Normally, availability is not specified alone. Both reliability 
and maintainability are usually specified together with availability. 
For durable, continuous operated hardware, reliability can be specified 
by mean-time-to-failure (MTTF). Given the exponential form of time 
to failure distribution, which is applicable to most types of complex 
equipment, MTTF is a constant, the reciprocal of the failure rate. 

Maintainability can be specified by mean-time-to-repair (MTTR) 
and mean-preventive-maintenance-time. Corrective maintenance 
times are often found to be well described by the exponential distribu- 
tion. Given that the time-to-repair is distributed exponential, MTTR 
is a constant, the reciprocal of the repair rate. Accordingly we 


make the following definitions. 


(1) A(t) - Availability at time t, or the probability that 
an item will be operable ata stated instant in 


time. tZ 0 


(2) IR - Interval availability is the time average of A(t) 
during intervals of length T, andis readily 
obtained from the equation for the average value 


of a function. 





- 
<a A(t)d > 

ies. (t )dt 

A, = (2) i 

i = a f A(t) dt 

J ° 
c 

Alternately, it may be considered that there is 

some probability distribution h(t) on demand time. 


Then the interval availability is given by 


oO 


Epa An (tle f aw h(t) dt 
oO 
In the special case where h(t) is uniform on the 


interval [0,T] 


h(t) = O42teT 


0 tL 
T 
= yl 
then A, - 7 J A(t) dt. 
O 
(3) A - the steady-state availability is the limiting interval 


availability as T —> scm 
vs 


A = lim xf A(t) dt 
O 


T > 
It can also be shown that 


A = lim A(t) 
TT co 


If the components do not undergo checkout or repair its interval avail- 
ability approaches zero. For systems with exponentially distributed 


failure times where 





Thus the steady state or long term uptime ratio of a component which 
does not undergo repair approaches zero. 

If the component is subject to repair, its steady state availability 
is not zero. Assume that both failure and repair times are exponentially 
distributed with means e and E: respectively. Given that A(0) = 


a 


i.e., component available at start of mission, then 


ue + jn 


repailLe rave 


A(t) 


i 


A = failure rate 
ee lent) = =o 
Vee IN 
t—> oO 


7 MTTF 
MTTF + MTTR 


Catron [2] has shown that exact values of A(t) can be and were 
computed for repair distributions other than exponential. The use 


of interval availability as t—> co i.e., 
7 


— 
A. = = fama>a 


4) 
to approximate the true availability over the first mission time as is 


currently being done in practice, is a conservative procedure. 
The amount by which the limiting interval availability under- 
estimates the true availability during the first mission time is 
dependent upon the limiting value, i.e., the lower limiting value, 
the more it underestimates the true availability. 
(4) Ag - System availability. For aseries system the 


system availability is defined as: 





A. = [I A. 
S : 1 
Ns 


~, 


where A. is the availability of the component of 


aye mea mc es °K 


Bb. PROBLEM 

The problem considered in this paper is to check the accuracy of 
a if 
Ag Tee or Ag based 


upon component failure time data and repair time data. 


a proposed lower confidence limit proceaure 


(1) When assuming exponential failure rate of. 
and exponential repair rate he Ene A. is given by 
fA 
oak = patio 
: [*- + Xe 
The following ranges of the failure rate ( XK.) and repair rate 


(9; has been used in the simulation. 


0.0054; <= 0.01 


3.04 7 


which implies that the ranges of MTTF of component i, 
denoted by K : and MTTR of component i, denoted od sia 


are as follows. 


1004 k » & 200 


oo! ea es 
y= 


(2) The equation of component availability can be writtenas 


yee = = i aa 
; [EF KE ei ower ee 
1 e 


but “A. = TF mepe = 1 a 1 pid = (xs es) vee 


8 


|e) 
Uo] — 





provided the series converges, which it does since Both loa 
Since the product ¢@ 44; is very small, in this case 


0.00054 <i 2 0.00333... 


Since oa ae will have small variance we shall apply the 


central limit theorem, and fit a normal distribution to 
KK 
Ro = Zz oy 
c=l 


Kain 
Then Ag = 1- 2 Xe pe . P 
cI, Ag 2 





Ni 
Assuming that Ag ey (7; eee 
° J 
oR. 
aD Ne 
Kn i cS ee axe 
7 Xv ¥ 
implies that Act gat Ag r (2. e 7 me 


C. 2a aRNAL IVE METHOD 
An alternative method was investigated, which fitted a 2 para- 


meter gamma distribution by method of moments to the distribution of 


Nn 
A 
A K eee : 
Ag= If he = 1 ee 
b= C=] (2: - 


/ /\ 
where x : sila wat are theveetimates of x : and f : respectively. 


Under the assumption that time to failure and time to repair of 
the components were distributed exponentially with failure rate % 


and repair rate /9 : respectively. The following approach was used 


to make the fit: 





barge Zt a) 


c=) 





sat Kit [Be 
Q.=1-A.=l- A 


dit (3< 


Thus by Taylor's expansion 


ke te 
a lsh J -lua-e, : Z. (On: 0% ) 
5 








Ve 
=Z Ye Be 


and 


.O> 
r >> 
X > 


os 


where  . and fP , are the estimate of & . and [s respectively 


i 
The 2 parameter gamma, [ ie Y) was fitted by method of 


Ys 
moments to the distribution of Ye 


TaN 
Bt Ys ) = wa 





A a 
a Yo 
var 7" 7 
me 4 
Vary) ae 
men, oy y.) | 
_ al 
Var (Ys) 
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by replacing = x. and -?; by be , and ia, 
LE reay 2 Y y. is distributed ue 
> ar 


From this assumption a one sided lower 100(1 - KX )% C.L. for Ag 


was computed. 
A 2 
Pek, P2V Ws ar) =1— co 
A r 
= i PE 
P(e ¥ 2 exp (“7s 1? 1\\ = | — of 
A 
P(AgZ Asi (x)= 1- 
A Oe 
where [2 }7~ ] = smallest integer greater than or equal to 2 I~ 


The procedure showed to be inaccurate due to small values of [~ 


SB PROCEDURE 
es 
The procedure used to fit a normal distribution to 1- Zi lec 
as described in B is based on computer simulation. 
The computations of the estimates of Xx a and f? is (lade 2° os 
are based on the assumption that we are givena series system with 


K types of components. 


a - time to failure of the j-th component of type i. 
The time to failure is assumed exponential with 
failure rate X rc 
ae time to repair of the j-th component of type i. The time 


to repair is assumed exponential with repair 7 rt 


VecatGel.. are independent. 
1) ij 


ell 





(1) A random sample of size N. components of type i are tested 


until failure or a specified planned test time (PTT). 


ee is the PTT of the j-th component of type i 


if if the j-th component of type i fail, i.e., 
life time less than fe 


~ 
oT 


“ O if the j-th component of type i lives beven 
ae 
N. 
Z, x oa) = # of failures of type i components. 


The components were tested against the (PTT) Left and the observed 


operating times ra set equal to ee or or whichever is the smaller. 


(2) A random sample of size n. of failed components of type i 


are repaired with repair times Ty T.5: sf, in, . 


(3) Using the test data described in C 1 and 2 the estimates of 


KM. ae were computed. 





Jc ee 
NM: a 
J=l 
Aw Be 
= € iq 
fee ‘eae 


(4) With the data in C-3 1000 values of Ag Lo ) were computed. 


A la K a je x 
oe Ake : 
A i ee > = wa 


Sie) ~ - wo 
=| Z Tei; x 
Bl ae 
ran lca. he 
where Ag =) are vA so aaa and 
C=/ 


2 





- is the 100(1- &) th percentile point of the N(0, 1) distribution. 


BE. ANALYSIS 


The generated values of A. L( ) was ordered from lowest to 
highest and the 1000(1- ec )% one was compared to Ag. 
The An Tee) was also compared to Ag and the probability of 
success was computed by the formula 
‘~ 
P(Asi(x ) = As) 
A 


based on the 1000 values of Ag L( of ) 


The results are shown in Appendix A for the different combinations 


” N,, n. and eget To avoid fluctuations due to the 


a 
of NX GB; 
{ a 


random number generator, the values of Ag L(X ) at different 's 


was computed by only changing Z x 
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MT, Sec CURACY RESULTS 


By changing the values of the input parameters the accuracy of the 
procedure was checked. The following conclusions were made based 
on the simulation: 

a. Use of the steady state availability of the system (Ag) to ap- 


A 
proximate the estimated availability A ) over the mission time 


SLOK 
is a conservative procedure. 


b. The amount by which the steady state availability differ from 


the estimated 100(1- & )th percentile point is very small. 


Kixample 
Caceclo=sppendix A. 
With the given parameters the systems steady state availability is 


computed, giving 


kK , 
aii a = .9925 
ix | [Bc + CGE 
A 


The computed 95th percentile point of the distribution of Ag OLS 





is equal to .9931 
A Z _ , 0006 
ines, Agee Nereis) | = | . 9925 = EI | ena 


Similarly the 90th and 80th percentile points of the distribution of 


A ix 
Ag et) and Ag Th) are .9930 and .9929 respectively. By increasing 


the number of items which are life tested from 50 to 100 and keeping 
the other parameters constant, we get the result given by case 5. 
That is, 


Ag = .9925 (as before) 
‘a = . 9928 
See. O5)0 a 
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This implies that 
A 





As - Agr 950)| * |. 9925 - .9928| = 
e /\ “A « . . 
In this case Asi 900) and Ag 800) > ° 9926 which implies that 
A 


SO0001 





Bg EES Thine )| as 


for (X = 0.1 and of = 0.2 


a 


The true levels of confidence associated with A as a lower 


> Coy) 


100(1- & ) confidence limit for A, are given in columns 10 through 12 


S 
/\ 
for X =.05, .10, .20 respectively. That is for Case 3, Ag L(. 05) 


is an 90. 8% lower confidence limit for A.. rather than 95% lower 


S 


confidence limit. This is a measure of the inaccuracy of Ag L(. 05) 


as a lower 95% confidence limit for A, for the given parameter values 


AN 
in Case 3. Likewise Ag L(. 10) 


limit and A, L(. 20) #5 really a 76.9% lower confidence limit. 


This variation between the true level of confidence and the proposed 


S 


is really a 84.5% lower confidence 


level of confidence (e.g., 90.8% vs. 95%) are due to the very small 
IW 


variance of A for the cases considered. It is felt that for 


SL( XK ) 
these cases the quantities 


I~ 


aS - Cie ~~ ) 


are better measures for the accuracy of this procedure. 
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1% 
Derivation of Mean and Variance of Ag 
thf “-# az 
TT ie = ‘ ar. 
ws : iz (e CK sf ae sear 


The MTTR should always be less than one mission unit, which implies 


Ki Kis By 


that 


thus 
| 
A. = “1 ee. = | 105. Oe = (2g. 2 + eceees 
1 Si Me fac * (pec) s 
provided the series converges which it does since Xe fee Gs ie 


Since the product Oe mes Very Simei 


Thus 


i te 

=~ s A —- er y Ke (i 
oa i \ ~dipi)= | Best J 

c= 


Since : and (>; will have small variance we shall apply the central 


limit theorem, and fit a normal distribution to 


Ke 
Ag = | - we oO & 
oS! 
Thus IN te Goa 
A ae /n 


when Lee is the estimated MTTR and &% ; is the estimated failure 


rate. Based on the assumption that time to failure and time to 
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repair is distributed exponentially we can find the mean and variance 


/\ 


of Ac. 


aN K \_ Saas _, “va a 4 A 
BIA.) =1-ZE(X)Elpi) = [—~ZAipé ee 
a 7 cel 
| n 7 
By using the equations given by OD 29304 [3], forcXy : and Var (& ) 


~~ 
we can find the variance for Ag 





Me) ve # of failures ZN 
— test time EN | 
m i 
Bee 3 )= XX. 
A ae 
Va~ (Xf .) = = 
i 
@ te 


= sum of all — test times on et a 
i (in mission units). 


= 

a5 

@ 

Ky 

@ 
i iS 
a 
a. 

| 





4 fe 
: = # of repair times — ee leg 
y at = 
U]- 


ke 


var( £2.) = ab 
Thus Var (Aq) = Var (1 - = a: i fe i) 
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